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An integrable Kondo lattice model, which describes a strongly correlated electron host interacting 
with a spin-1/2 lattice, is proposed. It is found that with the variations of the Kondo coupling J, 
the hole concentration and the magnetic field H, the system may fall into a variety of phases. 
The phase boundaries of the ground state are determined exactly. The marginal excitations and the 
quantum critical behavior at the phase boundaries are discussed. 
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Metallic compounds containing partially filled /-bands belong to the general category of heavy fermions Re- 
cently, with the discovery of the Kondo insulators |^ and the non-Fermi-liquid behavior, |^ the interest in this field 
has been greatly renewed. Especially, the non-Fermi-liquid behavior found in some heavy fermion compounds, which 
stimulates a strong challenge to Landau's Fermi liquid theory, reveals the quantum critical behavior of these sys- 
tems at low temperatures. The heavy fermion systems are usually modeled by the periodic Anderson model or 
the Kondo lattice model in some limiting cases. The single-impurity Kondo problem has been studied extensively 
and exact solutions |^,^ were obtained. In addition, important progress has been achieved recently for the Kondo 
problem in strongly correlated hosts [0-0. Nevertheless, the understanding to the Kondo lattice systems is rather 
unsatisfactory. Though many efforts have been made, an exactly solvable Kondo lattice model, which may provide 
us some crucial information of a heavy fermion system, is still absent. In fact, a generic two-impurity problem is very 
hard to be solved even in one dimension. I note a few integrable models consisting of many impurities have been 
proposed |ri|,|l^ . The impurities in these models are very artificial and are almost independent each other since only 
forward scattering between the conduction electrons and the impurities survive. The physical effect of these impurities 
is additive and therefore the problem is still at the single-impurity level. 

In this letter, I propose an exactly solvable Kondo lattice model consisting of a correlated electron host interacting 
with a spin-1/2 lattice. In one dimension, I show the model is exactly solvable via algebraic Bethe ansatz. The model 
hamiltonian reads: 

H = t'^his+5{i + n ■ n+s) + j'^St ■ fi, (1) 



-{n^ + Uj) + 1, 

where c|^^ (ci.o-) are the creation (annihilation) operators of the conduction electrons with spin a on site i; Si = 

\ So- ct' 4 CT^o-,o-'Ci,CT' are the spin operators of the conduction electrons (f the Pauli matrices); rii denotes the electron 
number on site i; V means the constriction Ui <\]Ti are the Pauli matrices indicating the local spin on site i; t and J 
are two real constants indicating the hoping amplitude and the Kondo coupling constant, respectively; i — (ii, - ■ ■ 
in d dimension and S are the basic vectors of the lattice, 5 = 1 in ID and (1,0), (0, 1) in 2D etc.. Obviously, Eq.(l) 
describes an S'f7(3)-invariant t — J system coupled with a spin lattice. The first term of Eq.(l) represents the 
hoping and interactions of the conduction electrons, which depend not only on the local electron states, but also on 
the local spin environment. As discussed in some previous works, p^ ] the electron- spin interactions contained in the 
first term can be mediated either by phonons or ordinary Coulomb repulsions. The second term describes the usual 
Kondo interactions. 

In the present form, the solvability of Eq.(l) is rather hidden. We note that the first term can be rewritten as 
g Pi^i+s, where P^.j is the permutation operator between the i-th site and the j-th site. For any given orthogonal 
and complete set of Dirac states {\ai >} spanning the Hilbert space of the i-th site, Pij can be expressed as 
(3 -^t^ -^j" ' where X"^ = \ai >< (3i\ represent the Hubbard operators. Now we check the possible states of a 
single site. A natural choice of {\ai >} is |7i,Tf >, where 7^ —1,0,1 denote one electron with spin up, no electron 
and one electron with spin down, respectively; and —'\, | denote the two components of the local spin. Obviously, 
the local Hilbert space is 6-dimensional and therefore the first term of Eq.(l) is iS'[/(6)-invariant. To diagonalize the 
whole hamiltonian, I introduce the following notations 
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|o>-^(| T,i>-U,T>), 

|1>=|0,T>, |2>=|0,i>, |3>=|T,T>, (2) 

|4>=^(| T,i>+U,T>), |5>=U,i>. 

|0 > represents a Kondo singlet; |1 >, |2 > denote two hole states and |3 >, |4 > and |5 > indicate the Kondo triplets. 
Obviously, < a\l3 >= 5a, p- Therefore, I can rewrite Eq.(l) as (up to an irrelevant constant) 

i,6 a,/3— i 

-jY^iXl'+Xf). (3) 

i 

Now it is clear that the Kondo interaction term represents one of the conserved quantities of the system. In fact, 
we have six types of conserved charges Na = X)i ^i"" which correspond to the number of local state |a > in the 
whole system. After this transformation, the hamiltonian (3) takes the exact form of an S'f7(6)-invariant spin chain 
model introduced by Sutherland ]l3t but with the Kondo coupling J as an effective field. In any dimensions, we have 
five branches of elementary excitations relative to a reference state \a >g— \ai > (X" • • • (X" \aN >. For example, if we 
choose |0 >g as the pseudo vacuum, we have two degenerate hole bands (corresponding to |0, t> and |0, |>) and three 
degenerate triplet bands. For the half-filled case, there is a critical value of the Kondo coupling J. When J > Jc, 
|0 >g becomes the ground state and the degenerate triplet bands are empty. Note no hole state is allowed in this 
case. Consider a single excitation upon the ground state with momentum k = (fci, • • • , kj_). The excitation energy is 
e{k) = 4:tJ2'^^^[cos{kd) — 1] + 2 J. Therefore, Jc = 'itd (t > is assumed). For J > Jc, the triplet excitations are 
massive while for J ^ Jc, these excitations are marginal and the system may show quantum critical behavior. For 
example, the low temperature specific heat and susceptibility at the critical point behave as C ~ T''/^, x ~ T'^/^"^. 
In fact, the massive spin excitations in the half filled case correspond to a Kondo insulator phase if we allow double 
occupation of electrons on a single site. This can be realized by replacing the S'C/(3) t — J term hi,j in Eq.(l) by the 
S'[/(4) or S't7(2|2) Hubbard term, jl^] In the latter case, the system exhibits also a charge gap which takes the same 
value of the spin gap A — 2( J — Jc) at half filling. For any non-half-filling case, we have two degenerate "hole seas" 
in the ground state configuration and the charge and spin excitations are always gapless. However, there is still a 
critical point Jc- When J > Jc, the triplet states are forbidden in the ground state. This critical value is generally 
filling-dependent and is hard to be derived in higher dimensions. The Kondo singlets behave as spin polarons. Because 
there is a finite gap to excite a triplet when J > Jc, we expect some type of condensation of the Kondo singlets at 
low temperatures. 

Bethe ansatz. In one dimension, the hamiltonian (3) has the same algebraic structure of that introduced by 
Sutherland [ p^ but with a very different physical interpretation. Not losing generality, we set 4t = 1 in the following 
text and choose the Kondo insulator state |0 >g as the pseudo vacuum state. The local states ja^ > (a > 1) can be 
treated as colored hard-core bosons with the single occupation condition X]q=o "^f " = 1- '^^^ Bethe ansatz equations 
(BAE's) for arbitrary symmetric or supersymmetric models have been given in Ref. [13]. In our case, they are 

/x(l) A/i >(1) ^(l) -Ah \(2) , I 

I 3 ~ 2 \ ^ TT ~ ^ TT ^3 ~ H +2 / . ^ 

(1)^1 11 ,(1) _ ,(1) ^, 11 ,(1) _ ,(2) _ W 



.^3 + t/ 1=1 ^ \ + * fe=l ~ 



k 



n ^3 ^ - » ^ TT TT -^J - \ - 2 

X(") \(") , ,■ 11 11 x(«) _ >(t) J ' 

1=1 ^fc + * t=n±l k=l ' 2 

n = 2,3,4,5; 



with the eigenvalue of the hamiltonian (3) as 



Ml 1 

^=-E(^TtI 2J)^JiN,+N2), (5) 



where M„ = Nn + ■ ■ ■ + and Mq = 0, A^"'' are the rapidities of the flavor waves (holons and spinous). Notice that 
the boundary condition X"^ = AT^'^^ has been used in deriving Eq.(4) and an irrelevant constant has been omitted 
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in Eq.(5). Relative to the Kondo insulator state |0 >g, the unoccupied states |0,t> and |0,J,> can be treated as 
Kondo holes which are responsible to the dynamical properties of a doped Kondo insulator. For large enough J > J^, 
there is no triplet states in the ground state configuration, i.e., M3 = 0. The effective low-energy hamiltonian of the 
system is therefore equivalent to an SU (3)-invariant t — J model. 

The phase boundary Jc{nh)- For a given Kondo-hole concentration = {Ni + N2)/N, there is a phase boundary 
Jc{nh) above which the spin triplet states will be eliminated from the ground state. Set Pn(A) as the distribution of 
A^"^-modes in the ground state. When J > Jc, from Eq.(4) we have 

/A2 pAi 
ai{X- v)p2{v)dv - / a2{X- i')pi{v)dv, 
-A2 J — Ai 

/Ai /'A2 
ai(A - v)pi{v)dv - / 02 (A - v)p2{i')dv, (6) 
-Ai J —^2 

/On(A) = 0, n>3, 

where a„(A) = n/[27r(A^ + n^/4)] and the cutoffs Ai^2 are determined by 

r-Ai 



pi{X)dX = nh, / pi{X)dX = -nh. 

-Ai J-A2 ^ 



(7) 



By integrating the second equation of Eq.(6), we obtain A2 = 00. The excitation energy of a triplet mode can be 
exactly derived by considering the process Mi Mi + 1, M2 M2 + 1 and M3 = 1. Such an excitation can be 
realized by adding a A^^^ mode Ap^^ above the A^^^ A(3)-mode A^^^ and a A^^) hole A^ ^ in the A^^) -sea. After 

some manipulation we obtain the excitation energy /^E{x'p\x^^\x^^) as 

A^ = ei(AW)-e2(Af ) + e3(A(3^), (8) 
where the dressed energy en(A) satisfy (see, for example, ref.[16]) 



/oo 
ai(A - v)e2{y)du 
-00 



i: 



(■A, 

a2(A - i/)ei{i/)du, 

-Ai 

/•Ai />oo 

e2(A) = / ai{X — i')ei{v)di' — / a2{X — v)e2{v)dv, (9) 

-Ai J —00 

POO 



/oo 
ai(A - i')e2{i')di', 
-00 



where = — 7rai(Ai) denotes the chemical potential. The energy gap associated with this excitation is given by 
^{nh) = £3(0) (Notice that ei(±Ai) = and e2(±oo) = 0). A is a monotonically increasing function of Uh and ranges 
from 2(J — 2) for n;i = to 2J for Uh = 1, while Jc is a monotonically decreasing function of Uh and ranges from 
Jc = 2 for ri/i = to Jc = for n/j = 1. For Uh = 2/3, the dressed energy e3(A) reads: 

«(A,„. = 2/3) = --yj-j^<!. + 2X (10) 

The energy gap A (2/3) and the critical value Jc(2/3) can be easily derived as 

A(2/3) = 2J-^ + iln3, Jc(2/3) = ^ - J ln3. (11) 

Quantum critical behavior at J = Jc{nh). When J < Jc{nh), the triplet excitations become massless and the 
system behaves as a Luttinger liquid with a holon band and four spinon bands. Exactly at the critical point J = Jc(n/i), 
the triplet excitations are marginal, indicating a quantum phase transition at this point. To see it clearly, let us consider 
the dispersion relation of the triplet excitations for J = Jcijih)- From the third equation of the BAE's (4) we know 
that a single A^^^-mode with rapidity A is quantized as 
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3 = 



where / is an arbitrary integer or half integer depending on the parity of Ah- Therefore, the left hand side of Eq.(12) 
can be treated as the quasi-momenta fc(A) of the A'^^^-mode. [ p^ In the thermodynamic hmit N oo, 

/oc 
arctan[2(A - i^)]p2{i^)diy. (13) 
-oo 

From Eqs.(6) and (9) I find that the velocity of the A^'^-' mode — hm^^^o de3(X)/dk{X) = 0, implying a finite mass 
of this excitation. Therefore, the dispersion relation takes the form: e3[fc(A)] = fc^/(2m), for fc 0. The effective 
mass of the excitation reads 



m = Van 



(14) 



It can be easily demonstrated that m takes positive values for any given ri/j. At very low temperatures, only few 
excitations exist and behave as a quasi-ideal quantum gas obeying the Pauli exclusion principle but with a zero effective 
chemical potential. The quantum critical behavior at J = Jc{nh) is mainly governed by the marginal excitations. For 
example, the low-temperature specific heat and susceptibility of the system behave as 

C^T^ (15) 

The divergence of the susceptibility at T = is due to the singularity of the density of states of the marginal 
excitations. 

Response to the magnetic field. When J < Jc{nh), the system behaves as a five-component Luttinger liquid, 
while for J > Jc{nh), an external magnetic field may drive some quantum phase transitions at zero temperature. For 
a given hole concentration n/i, there are three critical fields H^, and H^. In a very weak field, the susceptibility is 
Pauli type due to the response of the Kondo holes. When H reaches H^, the Kondo holes are completely polarized, 
while when H = H^, the Kondo singlets begin to be polarized and the zero-temperature susceptibility has a singularity 
at this point, x{H) ~ {H — H^)^^ {H > H^). For a strong enough H > H^, both the Kondo holes and the spin 
singlets are completely polarized. ranges from (n/,, 0) to 4 {uh = 1), 2(J — 2) < < 2 J and ranges from 
2( J -I- 2) {uh = 0) to 4 {nh — 1). For < H < H^, a magnetization plateau occurs since in this case, the Kondo 
holes are completely polarized while H is still not strong enough to excite the triplet modes. For < H^, singlet 
polarization occurs before the saturation of the holes' magnetization. There is no magnetization plateau for H < 
but there are a singularity aX H = and a kink at H = in the x~ H curve. To see the situation clearly, let us 
consider the ground state properties of nh ~ 1/2 and J > Jc(l/2) case. With a magnetic field, the eigenenergy reads 
(up to an irrelevant constant) 

Ml 1 

Si 



^ = - E ^TtI 2 JiVo - F(7V3 - iVs) - \h{N^ - N2). (16) 
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In this case, < < H^. When H = Hi, the ground state configuration is described by A'^o = = A^'i = Mi 
and Nn = for n > 2. The density of A'^-* is still given by Eq.(6) but with /02(A) = and Ai = 00. The critical field 
Hi can be derived by considering the excitation process A^i N/2 ~ 1, N2 1. This is realized by putting a hole 
m the A(i) sea and adding a particle to the A^^^ band. Denoting the rapidities of the hole and the particle as A^ and 
Ap, respectively and setting 6pi{X)/N as the change of pi(A) due to A/j and Xp, from Eq.(6) we have 



^pi(A) = - / a2(A - v)5pi{v)dv 

J —00 

+ai{X-Xp)-6{X-Xh). (17) 
The excitation energy associated with this excitation is 

e{Xh,Xp) = -Tr ai{X)Spi{X)dX + H. (18) 
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Solving Eq.(17) by Fourier transformation and substituting it into Eq.(18), we readily obtain the energy gap Ai = 
e(oo,0) = TJ — ln2. Obviously, TJ^ = ln2. can be derived in a similar way. For convenience, we choose |1 >g 
as the vacuum state. The excitation breaking a Kondo singlet corresponds to Mi — Nq N/2 — 1, M2 1. The 
energy gap associated with this excitation is A2 = 2 J — iJ — In 2. Therefore, = 2J — In 2. To derive H^, we 
choose still |1 >g as the vacuum state. When H > H^, the ground state is described by A^i = A^3 = N/2. With the 
same procedure we readily obtain the energy gap associated with this excitation reads A3 = — 2 J — In 2. Therefore 
Hi = 2J + ln2. 

Thermal BAE's. The thermal BAE's can be derived by following the standard method. | p7| , p^ Here I give the 
result and omit the details. After some manipulation, we obtain the density of the free energy as 



/ = -^E / fln(A)ln[l + ,7-;(A)]dA, (19) 



=1 



where rin,aW are the elements of the following integral equations 

G[ln(l + ?7„+i,„) + (1 - ln(l + ?7„_i^a)] 

= ^'5„,i<S„,i + G[(l - 5„,5) ln(l + 77-^+1) 



with the boundary conditions 



-(1 - <5„.i) ln(l + r,-i„_i)] + ln7M,„, (20) 
n = l,2---; a = l,2---5. 



limi^^^Ca, (21) 



where G is an integral operator with the kernel g{X) = 1/2 cosh(7rA); (i = J/T - n/T - H/2T, C2 = C4 = Cs = H/T, 
(■3 — J/T + /i/T — 3i//2 and /i denotes the chemical potential of the Kondo holes. Eqs. (20-21) can be solved with 
asymptotic expansions |^,^ or numerically. Details will be given elsewhere. 
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